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Summary 


A program is underway to improve the efficiency of a three-dimensional, unsteady, Navier- 
Stokes code and generalize it for nozzle and turbopump geometries. Code modifications will 
include the implementation of parallel processing software, incorporating new physical mod- 
els and generalizing the multi-block capability to allow the simultaneous simulation of nozzles 
and turbopumps. The modified code will be applied to relevant nozzle and turbopump con- 
figurations. 

The following results are described in the present progress report: 

• Details of code modifications 

• Numerical results for a nozzle configuration 

• Numerical results for interacting nozzle/blade/vane geometries 

• An update on the parallelization of the Navier-Stokes code using MPI software 
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Introduction 


Flow unsteadiness is a major factor in turbine performance and durability. This is es- 
pecially true if the turbine is a high work design, compact, transonic, supersonic, counter 
rotating, or uses a dense drive gas. The vast majority of modern rocket turbine designs 
fall into these categories. For example, the (Space Transportation Main Engine) STME fuel 
turbine, a high work, transonic design, was found to have an unsteady interrow shock which 
reduced efficiency by 2 points and increased dynamic loading by 24 percent. The Revolu- 
tionary Reusable Technology Turbopump (RRTT), which uses full flow oxygen for its drive 
gas, was found to shed vortices with such energy as to raise serious blade durability concerns. 
In both cases, the sources of the problems were uncovered (before turbopump testing) with 
the application of validated, unsteady computational fluid dynamics (CFD) to the designs. 
In the case of the RRTT and the Alternate Turbopump Development (ATD) turbines, the 
unsteady CFD codes have been used not just to identify problems, but to guide designs 
which mitigate problems due to unsteadiness. Using unsteady flow analyses as a part of the 
design process will lead to turbine designs with higher performance (which effects tempera- 
ture and mass flow rate) and fewer dynamics problems. One drawback to three-dimensional 
unsteady flow analyses, however, is that computation times are often on the order of weeks 
or months. This time scale does not allow for parametric design studies. The objective 
of the current program is to generalize and improve the efficiency of a three-dimensional, 
unsteady, Navier-Stokes code. The improved code will be applied to nozzle and turbopump 
geometries. 


Code Improvements 

The following improvements have been made to the CORSAIR analysis during the last 
quarter: 

• The inlet boundary conditions have been modified to allow the modelling of discrete 
nozzles. This is accomplished by applying the boundary conditions on a point-wise 
basis. The inlet boundary conditions can now be specified to be solid wall, subsonic 
inflow or supersonic inflow. Presently, it is assumed that the nozzles are circular or 
elliptical in shape. It should be noted that the turbulence quantities are still being 
calculated based on the airfoil and endwall surfaces. Eventually, the solid wall portions 
of the inlet should also be used to influence the turbulence quantities. 

• The technique for specifying the initial flow field in nozzle simulations has been mod- 
ified. For normal blade row calculations flow quantities are linearly interpolated be- 
tween the inlet and exit. In nozzle simulations (where it is assumed that the throat 
is choked) the flow quantities are interpolated between the inlet and the throat, and 
again between the throat and the exit. Note that for choked flow the throat quantities 
are known, and the area ratio can be used to check the Mach number at the inlet. 
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Application of this technique has been shown to reduce the CPU time by more than 
50% compared to initializing the flow field by interpolating the flow quantities between 
the inlet and exit. 

• A low-memory version of the code has been written in which many of the subroutines 
not necessary for nozzle simulations have been eliminated. This, in turn, has allowed 
the elimination of considerable branching logic. For a given problem size the low- 
memory version of the requires 40% less memory, and is 30% faster, than the standard 
analysis for nozzle simulations 

• PLOT3D files can now be specified to be formatted or unformatted. 

• A comprehensive and updated version of the CORSAIR User’s Manual was completed. 


Numerical Simulations 


FASTRAC Nozzle Operating in RPl 


The relevant turbine nozzle geometry in this investigation is from the FASTRAC pro- 
gram [1]. Two views of the nozzle grid (constant i- and constant k-planes, respectively) 
are shown in Figs. 1 and 2. The computational grid contains 131 points in the streamwise 
(?) direction, 81 grid points in circumferential (j) direction and 21 points in the radial (k) 
direction. Thus, the computational grid contains a total of 222,831 grid points. The noz- 
zle is designed to have an inlet Mach number of M\ ps 0.23 and an exit Mach number 
of M 2 « 2.13. The inlet static temperature is 1595° R and the ratio of specific heats is 
7 = 1.108. The Reynolds number was set at Re = lx 10 6 . 

Mach number contours at the inlet and exit of the nozzle are shown in Figs. 3 and 4, 
respectively. The average Mach number in the inlet plane is approximately Mi = 0.23, 
which is in excellent agreement with the design intent. The average Mach number in the 
exit plane is approximately M 2 = 2.18, which compares favorably with the design intent. 
Figures 5 thru 9 illustrate Mach contours along radial sections of the nozzle, beginning near 
the centerline and progressing towards the outer casing. These figures indicate that the 
highest Mach numbers are located near the centerline. Closer to the outer casing (Figs. 8 
and 9) the Mach number decreases due to the endwall boundary layer. Figure 10, which 
shows Mach contours in a circumferential ( j = 20) plane, highlights the thickening of the 
endwall boundary layer downstream of the throat. 

The predicted nozzle loss, which is defined as 


u = 


Pt inlet P^exit 


Pt 


inlet 


(i) 


was u> = 0.048 using the third-order accurate upwind scheme [2]. It is interesting to note that 
repeating the calculation using a second-order accurate upwind scheme the predicted loss 
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was u> = 0.061. Although the losses differed with order of accuracy, the average inlet/exit 
Mach numbers and flow patterns were similar. 

FASTRAC Nozzle Operating in Nitrogen 

A second simulation, which is nearly converged, was performed using nitrogen instead 
of RP as the operating fluid. Operating in nitrogen, the nozzle is designed to have an inlet 
Mach number of Mi « 0.23 and an exit Mach number of M 2 & 2.40. The inlet static 
temperature is 523° R and the ratio of specific heats is 7 = 1.467. Again, the Reynolds 
number was set at Re = 1 x 10 6 . The grid topology is identical to that used in the RP 
simulation. 

Mach number contours at the inlet and exit of the nozzle are shown in Figs. 11 and 
12, respectively. The average Mach number in the inlet plane is approximately M\ = 0.22, 
which agrees closely with the design intent. The average Mach number in the exit plane is 
approximately M 2 = 2.42, which compares favorably with the design intent. Figures 13 thru 
17 illustrate Mach contours along radial sections of the nozzle, beginning near the centerline 
and progressing towards the outer casing. Qualitatively, the results are similar to those 
predicted in the RP simulation. Figure 18, which shows Mach contours in a circumferential 
(j = 20) plane, highlights the endwall boundary layer development downstream of the throat. 

The predicted loss using the second-order accurate scheme is u = 0.071, which is higher 
than in the RP simulation. The calculation is currently being repeated using the thrid-order 
accurate scheme. 


FASTRAC Nozzle/Blade/Vane Geometry 

Numerical simulations have been performed for the FASTRAC turbine configuration 
consisting of a modelled nozzle exit, a rotor row and a vane row. The actual configuration 
consists of 26 nozzles, 147 rotor airfoils and 75 vane airfoils. In the current effort, two blade 
count approximations have been tested: 

1. l-nozzle/6-rotors/3-vanes - Using this approximation it is assumed that there are 150 
rotor airfoils, and the rotors are scaled by 147/150 to maintain the pitch-to-chord ratio. 
The modelling of the nozzle is then exact because 6 rotor passages will completely 
encompass one nozzle. 

2. l-nozzle/2-rotors/l-vane - This scaling also assumes there are 150 rotor airfoils, how- 
ever the nozzles now are approximated as a continuous slit around the annulus (see 
Fig. 19). 

The grid densities (number of passages xi xj xk) for the two simulations are presented in 
Tables 1 and 2, respectively. The total number of grid points in the l-nozzle/6-rotor/3-vane 
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Grid Type 

Inlet (Nozzle) 

Rotor 

Vane 

O 

- 

6x91x21x31 

3x91x21x31 

H 

1x14x51x31 

6x64x21x31 

3x82x31x31 

Total Points 

22,134 

605,430 

414,129 


Table 1: Grid dimensions for the l-nozzle/6-rotor/3-vane simulation. 


Grid Type 

Inlet (Nozzle) 

Rotor 

Vane 

0 

- 

2x91x21x31 

1x91x21x31 

H 

1x14x51x31 

2x64x21x31 

1x82x31x31 

Total Points 

22,134 

201,810 

138,043 


Table 2: Grid dimensions for the l-nozzle/2-rotor/l-vane simulation. 


in the l-nozzle/2-rotor/l-vane simulation. The average value of j/ + , the non-dimensional 
distance of the first grid line above the surface was approximately 1.5 for the airfoils surfaces 
and 5.0 for the endwall surfaces. A value of y + ~ 5.0 indicates that increased spanwise grid 
resolution may be needed. 

The turbine is designed to have an inlet Mach number of M\ & 2.11. The inlet static 
temperature is 1291° R and the ratio of specific heats is 7 = 1.108. The Reynolds number 
was set at Re = 1.4 x 10 6 . The results described below are for the l-nozzle/2-blade/l-vane 
simulation, which has completed 7 global cycles. A global cycle refers to two rotor airfoils 
passing one vane airfoil. 

Fifty snapshots of the turbine flow field were recorded during the course of one global 
cycle for the purpose of making an animation. Figures 22 thru 26 show instantaneous Mach 
number contours at 0%, 20%, 40%, 60% and 80% of a global cycle, respectively. Figures 27 
to 31 illustrate instantaneous static pressure contours at the same 5 time increments, while 
Figs. 32 thru 36 show instantaneous entropy contours. Note that the small breaks in the 
contour lines at the zonal boundaries are due to limitations in the graphics package and not 
to interpolation errors in the flow analysis. The graphics package included in the flow analysis 
first plots the contour lines on the H-grid, then “blanks” out the appropriate portion and 
plots the O-grid contours. Thus, the H- and O-grid contours are generated independently. 
Based on the snapshots of the flow field shown in Figs. 22 to 36 the following observations 
have been made: 

• The absence of strong shock or expansion waves far upstream of the rotor suggest that 
the rotor is operating at the unique incidence condition. 


10 





• There is a detached bow shock which stands approximately 3-4% of the rotor axial 
chord upstream of the rotor leading edge. The shock extends into the suction side of the 
rotor passage where it interacts with airfoil the boundary layer. The shock/boundary 
layer interaction induces flow separation at approximately 15% of the axial chord on 
the rotor suction surface. 

• The separated flow region on the rotor suction surface extends axially to the trailing 
edge of the airfoil and circumferentially to almost mid-passage (also see Fig. 37). The 
flow within the separated region is subsonic, but the flow in the remainder of the 
passage is supersonic. The entropy contours indicate there is a significant amount of 
loss generation associated with the separated flow region. 

• There is minimal unsteady interaction between the rotor and vane airfoils because of 
the relatively large axial gap between them. 

• The flow’ along the suction surface of the vane stays attached until immediately up- 
stream of the trailing edge. There is a shallow separation bubble originating from the 
leading edge of the pressure surface and extending to approximately 30% of the axial 
chord (also see Fig. 38). 

• There is a normal shock at approximately 53% of the axial chord on the suction surface 
of the vane. The shock extends a maximum of about 25% of the distance across the 
vane passage; the flow does not choke at any time during the course of a global cycle. 

In an effort to understand the development of the losses within turbine, instantaneous 
total pressure contours were plotted in axial planes 30% of the rotor axial chord upstream 
of the rotor passage, midway between the rotor and vane passages, and 47% of the vane 
axial chord downstream of the vane (see Figs. 39 to 41). Upstream of the rotor the contours 
indicate the presence of relatively thick endwall boundary layers (see Fig. 39). The reason 
for the thickness of the boundary layers is probably two- fold: 1) the solid wall portions of 
the nozzle at the inlet are located adjacent to the endwalls, and 2) the spanwise grid density 
needs to be increased. Note that the free stream value of the non-dimensional total pressure 
is Pt/Pco — 9.12 (see Table 3). Moving downstream to the region between the rotor and 
vane it is observed that the total pressure is only about 10-20% of its upstream value (see 
Fig. 40). There are two regions of low total pressure associated with the rotor wakes, but 
also tw T o regions near the shroud endwall where the total pressure is significantly higher. 
The source of these high total pressure regions is currently under investigation. Downstream 
of the vane there are loss regions associated with the vane wake and the endwall boundary 
layers (see Fig. 41). 

Table 3 contains the circumferentially and radially averaged flow variables at the inlet 
and exit of each blade row. Note that at the inlet only the flow coming through the nozzle, 
and not the solid wall values, are included in the average. 
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Variable 

INLET 

ROTOR 

VANE 

M in 

2.11 

1.96 

0.47 

Min 

2.11 

1.42 

0.47 

Mout 

1.96 

0.47 

0.39 

Min 

1.96 

0.74 

0.39 

Tt 

1601° R 

1595° R 

1439° R 

ft 

x L m 

1601° R 

1466° R 

1439° R 

Ti 

± L out 

1595° R 

1439° R 

1437° R 

Ti 

J l out 

1595° R 

1464° R 

1437° R 

(Pt/Poo)i n 

9.12 

7.66 

1.53 

(«/ft 

9.12 

3.09 

1.53 

(Ptlp~)„„ 

7.66 

1.53 

1.40 

(P‘/P~) M 

7.66 

1.85 

1.40 


72.0° 

73.4° 

-6.3° 

ft in 

72.0° 

66.8° 

-6.3° 

®out 

73.4° 

-6.3° 

o 

1 

ftout 

73.4° 

—52.0° 

-4.6° 

4* in 

0.0° 

-7.5° 

-1.6° 

V 1 in 

0.0° 

-4.4° 

-1.6° 

4>out 

-7.5° 

-1.6° 

-0.3° 

4^out 

-7.5° 

1.3° 

-0.3° 


Table 3: Average flow quantities. 
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Implementation of Parallel Processing Software 


The MPICH version of the Message Passing Interface (MPI) software has been down- 
loaded and installed on a Silicon Graphics Origin 200 workstation with two R1000 processors 
and 320 MB of memory [3]. Sample programs were used to verify the implementation. Tuto- 
rials on the MPI protocols have been obtained [3, 4], and the necessary routine calls will soon 
be added to the CORSAIR analysis. It is anticipated that the simulation will be distributed 
to processors on a per passage basis. 


Future Effort 

Work during the upcoming quarter will focus on: 

• Continuing the implementation of parallel processing software into the computational 
analysis 

• Continuing the nozzle/blade/vane simulations 
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Figure 4: Mach contours in the exit (i=131) plane of the FASTRAC turbine nozzle - RP. 
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Figure 5: Mach contours in the k=4 plane of the FASTRAC turbine nozzle - RP. 
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Figure 10: Mach contours in the j =20 plane of the FASTRAC turbine nozzle - RP. 
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Figure 12: Mach contours in the exit (i=131) plane of the FASTRAC turbine nozzle - N 2 . 


25 













0.50 


- 0.50 

Y(in) 


- 1 . 50 H 


- 2.50 


- 1.50 


INLET 


EXIT 




- 0.50 


X(in) 


0.50 


jffp 

mM 

* 7 “ 


1.40 


1.24 


1.09 


0.93 


0.78 


0.62 


0 . 4 “ 


0.31 


0.16 


0.00 


1.50 


Figure 17: Mach contours in the k=20 plane of the FASTRAC turbine nozzle - N 2 . 
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Figure 18: Mach contours in the j =20 plane of the FASTRAC turbine nozzle - N 2 . 
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Figure 19: Inlet nozzle definition. 
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Figure 20: Grid topology for FASTRAC turbine (x-y). 
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Figure 21: Grid topology for FASTRAC turbine (z-y). 
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Figure 22: Mach contours at midspan - 0% cycle. 
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Figure 24: Mach contours at midspan - 40% cycle. 
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Figure 26: Mach contours at midspan - 80% cycle. 
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Figure 28: Static pressure contours at midspan - 20% cycle. 
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Figure 30: Static pressure contours at midspan - 60% cycle. 
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Figure 34: Entropy contours at midspan - 40% cycle. 
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Figure 37: Relative- frame velocity vectors in the rotor - 80% cycle. 
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Figure 38: Relative-frame velocity vectors in the vane - 80% cycle. 
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